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Abstract 

A chain-like or a ID crystal-like universe made of intersecting mem- 
branes in extra dimensions in a Randall-Sundrum scenario is studied. A 
background gravitational metric satisfying the right boundary conditions 
is considered to study the localization of the scalar, gauge and Kalb- 
Ramond field. It is found that the wave function for the fields are Bloch 
waves. The mass modes equations are calculated allowing us to show the 
zero-gap mass behavior and the mass dispersion relation for each field. 
Finally we generalize all these results and consider q— forms in the crystal 
membrane universe. We show that, depending to the dimension D, the 
form q and the dilaton coupling A the mass spectrum can be the same for 
the different bosonic field studied. 
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1 Introduction 



When studying aspects related to localization of fields in membranes, we nec- 
essarily need to understand the solutions of differential equations with eigen- 
values. This eventually leads to the study a Schrodinger like equation open- 
ning the possibility for interpretation and comparison between methods and 
aspects of physics of extra dimensions and those of condensed matter physics. 
In some sense, it is another way to link these two areas, despite the celebrated 
AdS/CDM duality. Of course, there is no duality linking two theories with 
different regimens of coupling constant, but the basic mathematics of particles 
subject to potential barriers is just the same as those we attack in our training 
through the several famous quantum mechanical books. The question is how 
far can we go in this comparison. Besides the most basic example, the process 
of localization of fields (the zero modes or bound states), there is the study 
of resonances where the massive spectrum of particles interact with the mem- 
branes through the potential they generate (l||T4] . This kind of model provides 
a possible solution to the hierarchy problem and show how gravity is trapped 
to a membrane and also has important phenomenological implications to the 



standard model and cosmology 15 - 17 . Another very interesting step in that 



direction is the model regarding a universe with extra dimensions in such a way 



that the membranes itself generate a kind of chain or ID-crystal 18 -20 . In 
such situation we can think about membranes in every direction in the bulk, 
with intersections. These are configurations already discussed in the context of 
string theory phenomenology, where we can generate families of particles of the 



standard model plus some corrections [21-25 . Moreover, in string theory and 



supergravity the antisymmetric tensor fields arise naturally [26H28 and play an 



important role in dualization 29 30 . In particular, they appear in the R — R 
sector of each of the type II string theories. These tensor fields couple natu- 
rally to higher-dimensional extended objects, the so-called D— branes, and are 
important to make them stable. Besides this, they are related to the linking 



number of higher dimensional knots 31 and to the mechanism of topological 



mass generation |32[|33|. The rank of these antisymmetric tensors is defined by 



the dimension of the manifold 34 . Beyond that, these kind of fields also play 
an important role in the solution of the moduli stabilization problem of string 
theory |35}|37| . Because of these aspects, it is important to study higher rank 
tensor fields in membrane backgrounds. In this direction antisymmetric tensor 
fields have already been considered in models of extra dimension. 

In general the q— forms of highest rank do not have physical relevance because 
when the rank increases the number of gauge freedom also increases. However, 
that can be used to cancel the dynamics of the field in the brane [38], leading 
to the study of the mass spectrum of the two and three-form, for example, in 
refs. [39] and [40] in a context of five dimensions with codimension one. The 
coupling between the two and three-forms with the dilaton was studied, in 



different contexts, in 41-43 . In these scenarios, some facts about localization 
of fields are known: the scalar field (0— form) is localizable, but the vector gauge 
fields (1— form), the Kalb-Ramond field(2— form) and the three-form field are 
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not. That happens because in four dimensions the vector field is conformal 
and all information coming from warp factors drops out, necessarily rendering a 
non-normalizablc four dimensional effective action. Beyond that, Kehagias and 
Tamvakis [44] showed that the coupling between the dilaton and the vector gauge 
field produces localization of the vector field, and Alencar et.ai, considered these 



coupling with a three-form field where a condition for localization was found 45 



More recently, the problem has been generalized to g-forms 46 47 



The above mentioned facts have influence on the phenomenology of the four 
dimensional world we live [48]. The point we would like to address here is that, 
in the case of membranes distributed along extra dimensions which are slices 
of AdS spaces, we can find Bloch like solutions to the equations of motion. 
Therefore, one can argue that these models can present a similar behavior of 
quasi-particles like phonons, magnons, or electrons in the presence of a periodic 
potential. The goal here is to extend previous works in this subject to study the 
behavior of bosonic form fields in the crystal manyfold universe. We go right to 



the point of analysis of the fields subject to the background described in 18 
With that in mind, the organization of the paper is as follows. In the first 
section we shortly review the background of the crystal manyfold universe. In 
the second section we make a preliminary discussion of the general Schrodinger 
like equations we must treat. In the third, fourth, fiftieth and sixtieth sections we 
discuss the behavior of the scalar field, the gauge vector field, the Kalb-Ramond 
field and general g-forms, respectively. We discuss for each case the Bloch like 
behavior of the wave functions and give numerical results for mass bands. In 
the last section we present conclusions, more discussions and perspectives for 
future works. 



2 The Crystal Manyfold Background 



In this section we give a short review of the background metric obtained in the 
crystal manyfold. Wc only consider the simplest case with one extra dimension. 
Following ref. [18], the conformal metric is defined as 



with 



ds\ = n 2 (z)(T lfll ,dx'"dx u + dz 2 ), 
fT 1 = KS{z) + 1 



(1) 



(2) 



where in the above equation K = (L) 1 [L is the AdS radius). The function 
S(z) satisfy the equations 

*M _,£(-!>*(,- JO, 



dS(z) 
dz 



= 1, 
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and are given by (I is the spacing between the branes) 



= . 2pl - z, for (2p -l)l<z< 2pl; 

* z - 2pl, for 2pl < z < (2p + 1)1; { ' 



In order to analyze the dilaton with a conformal metric ([TJ), we will use the same 



procedure as in ref. 12 where the metric with dilaton is 

ds 2 = e^^^dx^dx" + e 2B{v) dy 2 , (4) 
and after solving Einstein's equation we get the relations [44]. 



B(y) = A(y)/4; tt = -V3M*A(y). (5) 

We introduce the parameter b such that B(y) = (1 — b)A(y), where b = 3/4 
is the case with dilaton and 6=1 otherwise. To obtain the conformal metric 
(!]), we use the transformation dy/dz = e A (v)-B( y ) = e bA( y ) = in Q We 
obtain A(z) — \nil(z)/b, where A(y) — A(z). 



3 The General Case 

Here we focus on the prototype Schrodinger-like equation that will appear 
through all this work. We consider the solution with some details since it will 
be important for all cases considered. As mentioned before we can relate the 
potentials with the ones obtained previously in [12] . In all the cases the equation 
of motion has the following general form 

Jf +P ' {v) iy + V[v) ^(y) = m2 Q(yW(y)> ( 6 ) 

where P{y) — ^A(y), Q(y) = e~ 2bA ( y * > and V(y) = for all fields except grav- 
ity, in which V(y) = 2A"(y) — 2(1 + b)A'(y) 2 . We can transform this into a 
Schrodinger-like equation through the transformations 

? = /(»). W = eWW, (7) 

ay 

with 

f(y) = VW), ®(y) = eMP(y)/2)Q(yr 1/4 - (8) 

As stressed before, when we have an analytical expression for dy/dz the potential 
is given by 

e(z) G{z) f(z) 
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where f(y) = f(z). With this we arrive to a Schrodinger like equation 



dz 2 



U 



■0(z) = m 2 ip(z), 



where the potential U(z) is given by 

U(z) = cA"(z) + c 2 [A'(z)} 2 . 
Now, considering A(z) = lnSl(z)/b. We get 

(ft- 1 )' 2 c(fi- 1 )" 



U(z) 



(n- 1 ) 2 b (n- 1 ) 



(10) 



(ii) 



(12) 



The potential has been expressed in terms of f2 _1 due to the form of the con- 
formal factor. It is interesting to note that the form of the potential coefficients 
simplifies the solution. Taking f2 given before we get the final Schrodingcr-like 
equation 



m 



1 

b ' b 2 ) {S{z)+L f 



)C E.(-ip<H*-jQ 

J & (S{z)+L) 



ijj = 0. (13) 



This equation is very similar to the Schrodinger equation for a particle propa- 
gating in a one-dimensional periodic potential (ID-crystal). In order to solve 
it we have to consider two adjacent elementary cells in the region < z < Al. 
Focusing on the region < z < 21 we need to find the appropriate boundary 
conditions. First of all we must have ip(l+) = ip{l-) and ip(2l + ) = i/j(21_). 
Where we use the definition l± = lim^o^ i e). To obtain the correct boundary 
condition for the derivative we simply integrate the above equation from Z_(2Z_) 
to (2Z_|_) to obtain 



2c 
bjl + L) 
2c 



^(2z+) - v , (2z_) = -^m)- 

The equation for the first elementary cell is 



(%+ ^4> for 



{21- z 



0<z <l 

■f^yrlp for I < Z < 21 



(14) 



(15) 



and the above equations can be expressed in terms of a Bessel equation. For 
this we just perform the transformations 



{y/u^(u), u = m(z + L), for < z < I; 
y/v$(v), v = m(2l — z + L), for l<z<2l, 



(16) 
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to get 

m 2 *" + u^' + (u 2 -v 2 )^> = 0, (17) 

where v 2 = (| + |) 2 and the same equation is valid for the variables u and v. It 
is important to mention that the quantity v 2 is related to the kind of field we 
consider as we will show later. Therefore the solution of the above equation is 
given by 

_ fV«(Aff+(u)+Bfl--(u)), 0<z<l, u = m(z + L); 

\yfu(CH+(v) + DH-(v)), l<z<2l, v = m(2l-z + L). 

Now using the transfer matrix technique we can relate the constants C, D 
with A, B. In order to impose boundary conditions for the wave function and 
for its first derivative we define 

E = ^/uH+{u),E{l) = y/m(l + L)h+; 

F = VuH-(u),F(l) = y/m(l + L)hZ, (19) 

and 

G = V^H+(v), G(l)=E(l); 

I = VvH-(v), I(l) = F(l), (20) 

with the following properties 

j- z E(u)\ z=l = mE'(l) = -mG'(l); (21) 

j- z F{u)\ z=l = mF'(l) = -mI'(l), 

where the prime means a derivative with respect to the argument. After a bit 
of algebra we obtain 

V(u) \ ( E F \(A^ 

^£tp(u)J \mE' mF'J \B y 

Itlyj = {-mC -Li){d)- (22) 
With this the boundary conditions can be accounted by imposing 



E(l) F(l) \ (A 

mE'{l) + jrffaE® mF'(l) + ^F(l)) { B 

E(l) F(l) \ (C 

-mE'(l) -mF'(l) J \D 



Inverting the matrix in the RHS 



(23) 



E(D F(l) y'_ 1 f-mF'(l) -F(l)\ 

-mE'(l) -mF'(l)J m 2 {i + L)(Kb+^ hth~_ x ) \+mE'(l) E(l) J 

(24) 
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wc obtain 



/ h~h+_ 1 +h+h~_ 1 

~~ hZh+_ 1 -hth~_ 1 
2h+h+ , 



-2h~h ly _ 1 
hZ h+_ 1 -hth~_ 
(h~h+_ 1 +h+h~_ 1 ) 



K 



y Kh+_ 1 -hth~_ 1 Kh+_ 1 -hth~_ 

Similarly, we can solve for the next cell (21 < z < Al) to find 

/ K h t-i+ h t h Z-i - 2 KK-i 



(25) 



V 



Kh+_^-hth~_ 
2h+h+_ 1 
h~h+ ,-hth- ' 



Kh+_ 1 -hth~_ 

(K h i-i+ h iK-i) 



-hth; 



= K 



(26) 



and therefore we have four constants A, B, Aand B to be determined. In order 
to get this we impose the following conditions 



V(2Z) = e 2i *V(0), 
i>{Al) = e 2lql ^(2l), 
^(21+) = ^(2/_), 

^(2/ + )-V'(2i_) = "^^(20- 

The first of the above conditions can be written as 



(27) 
(28) 
(29) 

(30) 



e 2iql (h+ K)(^=(h+ K)(^) = {hi K)K(fj, (31) 



where we have used H(2l) = H(0) = h, resulting in 

B = h+K 11 + h-K 21 -e 2i " l h+ A ^ 
h+K 12 + h„K 22 - e^ l K 

and identically for the second condition 

g = h+Kn + h-K 21 - e 2 ^ l h+ 
h+K 12 + KK 22 - e 2i " l fc 

From the third condition we get 

Ah+ + Bh- =e 2iql (Aht +Bh~ 
and using the last two relations we find 

A = e 2iql A. 



A = J A. 



(32) 



(33) 



(34) 



(35) 



The last two conditions are identical to the contour conditions used in z = I 
and z = 31, therefore from them we obtain 



K 



KK 



KK 



A 
fA 



(36) 
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where K is the same matrix as K with h — > h. Using now the Eqs. (32), (33 1 



and ( 35 1 we get 



2iql \ A / j 



fe 2 iq i J = KK (37) 
and finally we arrive to 

{j v n v -i +j' I /-in I/ )(j I/ n„_i + j v -in v ) - 3v-i3u{ju-iju + 3n v -in v ) 



cos(lq) 



^(jun v -i - ju-in v )(j v nu-i - iv~\n v ) 



^{jvny^x - 3v-inv){jvfi v -i - 3u-\n u ) 



(38) 



The above equation gives us a dispersion relation for the mass, where one can 
find, for example, the mass gap of the system considered. In order to obtain 
numerical results we have to discuss the magnitude of the parameters involved. 
We have that mc 2 = (fi-c)j where we can compare units of mass and length. 
For this purpose we use he ~ 200AfeV r x 10 _15 m. In the case where / > L, the 
corrections for the newtonian potential for particles localized to a 3-brane at 



the intersections of branes with negative tension is approximately given by 18 



W,- G „^( 1 + ^), (39, 

where a and b are constants of order unity. With that we can have gravity in the 
observational bounds, i.e., -p- < 1mm. For this case Mp L ~ Af^LH (H is the 
number of branes in the crystal and L is the AdS radius). The fundamental scale 
should satisfy M„ ~ • TeV. If I ~ eV^ 1 <C 1mm we find H ~ 10 16 branes and 
M* ~ lOOTeV (jr ~ 100). That is the result we put in our numerical analysis. 
In the next sections we present the results for each case of bosonic fields in 
detail. 



4 Scalar Field on the Crystal Manyfold 

In this section we consider the scalar field case. Before analyzing the massive 
modes we must look for the localization of the zero mode. The action for the 
scalar field without the dilaton coupling is given by 

S = J d 4 xdzv^GG MN d M ^>d N ^ (40) 

and the zero mode is obtained by considering $ = x( x )i which is equivalent to 
set m = 0. Considering the metric in the form of Eq. ([I]) we get for the effective 
action 

S = J dzn z {z) J d 4 x V ^d^ X d,X, (41) 

and according to the solution J fl 3 (z) < oo we have a well defined four dimen- 
sional action. It is clear that the scalar field is localized without the necessity of 
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the dilaton coupling. This coupling will not change the localizability of the field, 
but it must be considered because of the gauge field localization. Therefore, due 
to consistency it will be include in the analysis. 

To find the massive modes we should solve the equation of motion. Here we 
use the metric defined by Eq. Q. As we have found before, the potential is 
written in a general form 

U(z) = c i A"(z) + c 2 i [A'(z)} 2 . (42) 

The prime means a derivative with respect to z, and to obtain that result we 
have used dy/dz — e~ biA , where the label i — 1 and 2 stands for the case with 
and without the dilaton, respectively, and are given by 

3 ( 3 AV3M 3 ' 

h = -, d = 

b2 = I, c 2 = (43) 




and we get 



2AV3M 3 \ 2 



v 2 = 2. (44) 



5 The Gauge Field in Crystal Manyfold 

Now we consider the gauge field living in the crystal manyfold. The action for 
this field is given by 

Sx = / d 5 xV^G[Y MlM . 2 Y AhM % (45) 

where Ym 1 m 2 = ®[M 1 Xm 2 ] ^ s the field strength for the 1— form X. The study of 
the zero mode is made using Xm 1 = Xm 1 (x^) and, as said in the second section, 
we must consider the metric given in Eq. ([!]) for the case without the dilaton 
coupling. In this case we get the effective action 

Sx = J dzn J d 4 x[Y^ l2 Y^^}. (46) 

Like in the case with just one brane J dzQ = oo, and we have a ill defined action. 
However, we can consider the coupling of this field with the dilaton, with action 
given by 

S x = f d 5 xV^Ge- X7r [Y MlM2 Y AhA % (47) 
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where A is the dilaton coupling. With this and using the metric defined in Eq. 
Q the effective action becomes 



J dzrti- x ^ J d^Y^Y^}. (48) 



The integration must be performed in the crystal manyfold. However, as men- 
tioned previously, the finiteness of the z dependence is reduced to the finiteness 
of integral in one cell. With the expression for il we see that this is reached if 
A > -1/V3A/ 3 . This is the same condition as that obtained for the localization 
of gauge field in the model with just one brane. In order to analyze the massive 
modes, we must come back to the metric defined by Eq. Q. The equations of 
motion are given by 

d Ml (V^GG AhP G M ^Y PQ ) = 0. (49) 

We use a gauge freedom to fix X y = d^X^ = and by using the transfor- 
mations g| = e~ 3A / 4 and || = e~ A respectively, we reach the potentials for the 
Schrodinger-like equation for the cases with and without the dilaton as before 



U(z)= Ci A (z) + ci[A'(z)} 2 , (50) 



where now 



3 / 1 AV3M 3 \ 

h - 4' Cl= (j + ^— J' 

b2 = 1, C2 = \, (51) 



and we get 




(52) 



6 The Kalb-Ramond Field in Crystal Manyfold 

Now we turn our attention to the Kalb-Ramond field. The action for this field 
is defined by 

Sx= J d 5 xV^Y MlM2M3 Y M ^ M ^ (53) 

where Ym 1 m 2 m 3 = d[M 1 ^M 2 M 3 ] is the field strength for the 2— form X. Again, 
we must consider the zero mode given by Xm^m?, — Xm 1 m 2 {x' 1 ), with an effec- 
tive action in the conformal metric ([I]) given by 

S x = J dzil- 2 J d 4 x[Y^ l2 ^Y^^}, (54) 
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and, as in the gauge field case, we have a ill defined effective action since 
f dzfl~ 2 = oo. Now we turn to the action with the dilaton coupling defined 
by 

Sx= J d 5 xV^Ge-^Y MlM2M3 Y M ^ M ^ (55) 

where again A is the dilaton coupling. With this and using the metric in Eq. 
Q the effective action becomes 



Sx = J dzrt-*-^ J d*x[Y^ 2 Y^}. 



(56) 



With the expression for ft we see that the zero mode localization is reached 
if A > 7/4V3M 3 . A gain we find the same condition as that obtained for the 
localization of Kalb-Ramond field in the model with just one brane. Now we 
analyze the massive modes coming back to the metric defined by Eq. Q. The 
new equation of motion is: 



d M (V^GG MP G NQ G LR e- X7T H PQR ) = 0. 



(57) 



Here, we can use gauge freedom to fix X^ iy = d fll X fll ^ 2 = and by using the 



transformations ^ 

ay 



" 3A / 4 and 



<iy 



e respectively, we have the following 



expression for the constants 6^ and c 



ci = 



AV3A/ 3 ' 



h = 1, c 2 
and we obtain the solutions 



v 2 = 0. 



1 

2' 



(58) 




(59) 



7 The q— Form Field in Crystal Manyfold 

Now we can readily generalize our previous results to a q— form in a p— brane, 
where p = D — 2. In a recent paper the authors have considered this issue and 



we are going to use those results here 46 . Just as before we must consider the 



cases with and without the dilaton coupling. The action is given by 



->x 



d D xV^G{Y Ml ... Mq+1 Y M i- M ^} 



(60) 



and considering Xm\ m 2 m 3 = Xm 1 m 2 M 3 {x ij ') with the metric of Eq. ([I]) we arrive 
at the effective action 



S x - / <fefi D - 2 <* +1 > J d 4 x[Y^ 1 ... flg+1 Y^-^}. 



(61) 
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1| , 1 , 1 , |_ 

-1 -0.5 0.5 

ql/% 



Figure 1: The lowest mass modes for the case without dilaton and the scalar 
field v = 2 

It is evident why the gauge field is not localizable in five dimensions. From 
the above we get the condition q < (D — 3)/2, which for D = 5 give us q < I. 
We also see that in higher dimensions we can have localized form fields without 
the inclusion of the dilaton coupling. However, if we want to localize all the 
form fields we must consider the dilaton coupling with an action defined by 

Sx= J d D x^Ge-^[Y Ml ... Mg+1 Y Ah - M ^] (62) 

and considering again Xm 1 m 2 m 3 = XM 1 M 2 M 3 {x tl ) and Eq. Q we get the effec- 
tive action 

S x = J dz^v- 2 ^-^ J d 4 x[Y^..^ q+1 Y^-^}. (63) 

From the above expression we get the condition A > (8q — 4p + 3)/4v / 3AP for 
localization. We must note that all the previous conditions are enclosed in this 
one. This expression is the same obtained for the localization in the case of one 
brane. To consider the massive modes we must consider the metric defined by 
Eq. Q, and the equations of motion for this case is given by 

d M (V^9 MNl ■ ■ ■ 9 M * N ^Y Nl ... Nq+1 ) = 0. (64) 
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Using the gauge freedom to fix X fil ... fiq iy — d v X v ... )1 = and the transforma- 
tions ^| = e _3j4 / 4 and ^ = e~ A we arrive to the following expressions for the 
parameters bi and Cj; 



3 fa 3 

h = 4' Cl= -U + 8 

62 = 1, C2 = (f-<z), (65) 



with a = (8q — Ap — 3)/4 — AV3Af 3 . From these the solutions are given by 

2a\ 2 

^- ? Y. (66) 



8 Results and discussions 

The important information one can get from the calculations in previous sections 
is the allowed values of mass for each field. That can be obtained by solving 



numerically the equation (38 1 for the cases with and without dilaton. For the 
case without dilaton, we solve Eq. ( [38] ) for the scalar field using v = 2, for the 
gauge field v — 1, and the Kalb-Ramond fold v = 0. The allowed masses for 
the scalar field are shown in Fig. [T] where the lowest mass modes are plotted. 
As one can see, the lowest mode or the gap of mass has its maximum value 
at q = 0, and for small values of |g| there are only four allowed mass values. 
The interesting feature about the dispersion relation equation is that due to 
the property of the Bessel function j_i = —j\ and n_i = — ni, the gauge field 
and the Kalb-Ramond field have the same mass dispersion as shown in Fig. [2j 
The main characteristic of these fields is that around q = no mass is allowed. 
In fact for m > 0.5 the existence of mass particles is restricted to values of 
?w 0.5tt/Z. 

When considering the dilaton, depending on the coupling the order of the 



Bessel functions in the dispersion relation can be non-integer. For A = 1/ \J (3M 3 ) 
Figure [3] shows the mass dispersion for the q — 0, 1, and 2 forms in panels o, 
b, and c respectively. For the scalar field Fig. [3)2. shows a different behavior 
from the Fig. [T]the main effect is now the existence of more mass modes around 
q = 0. A similar effect happens for the gauge field, whereas in Fig. [l] there was 
no allowed mass for q = 0, one can see two possible values of mass in Fig. [3f). 
The dilaton coupling breaks the degeneracy of the mass modes for the gauge and 
Kalb-Ramond fields. Now, the two spectra are very different (see Fig. |3j:). The 
main effect of the coupling is to change the order of the Bessel functions. The 
behavior of the lowest mass modes against the Bessel function order is shown 
in Fig. § 



12 



-1 -0.5 0.5 

ql/% 

Figure 2: The mass modes for both the gauge [y = 1.0) and Kalb-Ramond 
{v = 0.0) fields 
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a 



b 



^2 




4 



2- 



-1 



-0.5 





ql/n 



0.5 



Figure 3: The mass spectrum for fields with a dilaton coupling A = (3M 3 ): 
a) for the scalar field (v — 3.17), b) for the gauge field (y — 1.84) and c) the 
Kalb-Ramond field [y = 0.5) 



In the Figure (4jz it is possible to distinguish two regimes: for v > 1 the 
mass values increases, while for v < 1 the masses have a parabolic behavior. 
For q = 0.0 Fig. |4f> shows that below v ps 5/3 no masses are allowed. These 
results leads to the fact that controlling the dilaton coupling one can generate 
or suppress mass modes. A similar effect can be obtained by playing around 
with the dimension D and the form of the field q. For example, in order to get 
v = 5/3 one can have a dilaton coupling A = 3(3M 3 ) -1 / 2 /4 and the condition 
q = (2 + D)/2, i.e., only an even dimension is allowed to give an integer q 
form: for D — 2 we must have q=2. When considering A = 7(3A/ 3 ) _1 / 2 /4 the 
condition is q — (3 + D)/2, in this case only an odd dimension can give an 
integer q. In this case for D = 1 we can have q = 2. 

Another interesting aspect of the mass dispersion relation is the lowest mass 
behavior against the spacing between the branes as shown in Figure [5j In 
there, for several values of v the mass mode decreases as m(l) oc l/l. For large 
separations the are not interfering with each other and the system behaves like 
a single brane and no gap appears for any any q-form or dilaton coupling. 
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Figure 4: The mass dispersion against the order of the Bessel function for: a) 
q = 0.5tt/;, and b) q = 0.0 

9 Conclusions 

In this paper we studied aspects of bosonic fields subject to a gravitational back- 
ground generated by a brane crystal system, where the simple case of branes 
distributed only along one extra spatial dimension was considered. The behav- 
ior of the scalar field, the gauge vector field, the Kalb-Ramond field and general 
g-forms was analyzed. The dispersion relation for the mass was presented for 
each field mentioned above. Without the dilaton we noted that the Gauge and 
the Kalb-Ramond fields have the same mass spectrum when D = 5. However 
for a lower dimension, for example D = 3, the scalar field and the gauge field 
will have the same spectrum. This is an interesting result, since different fields 
can have the same allowed values of mass for bosonic particles. Even with the 
introduction of the dilaton coupling, different forms with different couplings can 
have the same allowed values of mass. Therefore, the main conclusion is the fact 
different bosonic fields can produce particles with the same mass. Another im- 
portant aspect, is the analogy with the problem of an electron traveling through 
a periodic potential. In that case, one can calculate the allowed values of the 
electron's energy where the fields used here make the role of constraints in the 
potential of the ID-crystal. That opens the possibility to relate physical effects 
in condensed matter physics to the physics of extra dimensions. 
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A natural extension of this work is to study fermions subject to the same 
background. Besides, it is very important to avoid physical space-time singular- 
ities due to the delta function like branes. In this case, we could propose smooth 
configurations mimicking the crystal. Such configurations can be obtained as 
approximated multi-kink solutions of a specific scalar field configuration. These 
are discussions left for future researche. 



Acknowledgments 

We acknowledge the financial support provided by Fundagao Cearense de Apoio 
ao Desenvolvimento Cientffico e Tecnologico (FUNCAP), the Conselho Nacional 
de Desenvolvimento Cientffico e Tecnologico (CNPq) and FUNCAP/CNPq/PRONEX. 

This paper is dedicated to the memory of my wife Isabel Mara (R. R. 
Landim) . 



References 

[1] D. Baz eia, A. R. Gome s and L. Losano, Int. J. Mod. Phys. A 24, 1135 
(2009) [arXiv:0708.3530| [hep-th]]. 

[2] 
[3] 



W 
[5] 
[6] 
[7] 



[9 
[10 

[11 



R. C. Fonseca, F. A. Brito and L. Losano, arXiv:1211.0531 [hep-th]. 

A. E. R. Chum bes, J. M. Hoff d a Silva and M. B. Hott, Phys. Rev. D 85, 
085003 (2012) | arXiv:1108.3"82l"1 [hep-th]]. 

M. S. Cunha an d H. R. Christiansen, Phys. Rev. D 84, 085002 (2011) 
|arXiv:1109.3486| [hep-th]]. 

C. A. S. Almeida, M. M. Ferreira, Jr., A. R. Gomes and R. Casana, Phys. 
Rev. D 79, 125022 (2009) |arXiv:0901.3"543] [hep-th]]. 

Y. -X. Liu, J. Yan g, Z. -H. Zhao, C . -EFu and Y. -S. Duan, Phys. Rev. D 
80, 065019 (2009) [arXiv:0904.1785] [hep-th]]. 

Y. -X. Liu, H. -T. Li, Z. -H . Zhao, J. -X. Li and J. -R. Ren, JHEP 0910, 
091 (2009) [arXiv:0909.2312| [hep-th]]. 

W. T. Cruz, M. O. Tahi m and C. A. S. Almeida, Europhys. Lett. 88, 41001 
(2009) |arXiv:0912.i029| [hep-th]]. 

Z. -H. Zhao, Y . -X. Liu, H . -T. Li and Y. -Q. Wang, Phys. Rev. D 82, 
084030 (2010) [ arXiv:1004.2181| [hep-th]]. 

H. -T. Li, Y. -X. L iu, Z. -H. Zhao and H. Guo, Phys. Rev. D 83, 045006 
(2011) arXiv: 1006.4240 [hep-th]]. 

W. T. Cruz, A. R. Gom es and C. A. S. Almeida, Europhys. Lett. 96, 31001 
(2011) [arXiv:1110.3~i~04| [hep-th]]. 



16 



[12] G. Alencar, R. R. Landim, M. O. Tahim and R. N. C. Filho, 
|arXiv: 1207.3054] [hep-th]. 

[13] R. R. Landim, G. Alencar, M. O Tahim and R. N. Costa Filho, JHEP 
1202, 073 (2012) ljarXiv:1110.5855| [hep-th]]. 

[14] R. R. Landim, G. Alencar, M. O Tahim and R. N. Costa Filho, JHEP 
1108, 071 (2011) I jarXiv: 1105. 5573| [hep-th]]. 

[15] C. Csaki, J. Erlich and J. Terning, Phys. Rev. D 66, 064021 (2002) [hep- 
lph/0203034] . 

[16] C. Csa ki, M. Graesser, L . Randall and J. Terning, Phys. Rev. D 62, 045015 
(2000) [hep-ph/9911406] , 

[17] C. Csa ki, M. Graesser, C. F. Kolda and J. Terning, Phys. Lett. B 462, 34 
(1999) [hep-ph/9906513] , 

[18] N. Kaloper, Phys. Lett. B 474, 269 (2000) |hep-th/9912125] . 

[19] N. Kaloper, JHEP 0405, 061 (2004) [hep-th/0403208] . 

[20] I. Oda, Phys. Lett. B 480, 305 (2000) [hep-th/9908104l . 

[21] G. Aldazabal, S. Franco, L. E. Ibanez, R. Rabadan and A. M. Uranga, 
JHEP 0102, 047 (2001) [hep-ph/0011132] . 

[22] L. E. Ibanez, F Marchesano and R. Rabadan, JHEP 0111 (2001) 002 
|hep-th/0105155| . 

[23] R. Blumenhagen B. Kors, D. Lust and T. Ott, Nucl. Phys. B 616, 3 (2001) 
[hep-th/0107138] . 

[24] R. Blumenhagen, M. Cvetic, P. Langacker and G. Shiu, Ann. Rev. Nucl. 
Part. Sci. 55, 71 (2005) |hep-th/0502005] . 

[25] D. Lust, Class. Quant. Grav. 21, S1399 (2004) [hep-th/0401156] . 

[26] J. Polchinski, "String theory. Vol. 1: An introduction to the bosonic string," 
SPIRES entry Cambridge, UK: Univ. Pr. (1998) 402 p 

[27] J. Polchinski, "String theory. Vol. 2: Superstring theory and beyond," 
SPIRES entry Cambridge, UK: Univ. Pr. (1998) 531 p 

[28] P. Van Nieuwenhuizen, Phys. Rept. 68, 189 (1981). 

[29] A. Sma ilagic and E. Spallucci, Phys. Rev. D 61, 067701 (2000) |arXiv:hep- 
|th/9911089] , 

[30] A. Sm ailagic and E. Spallucci, Phys. Lett. B 489, 435 (2000) |arXiv:hep- 
|th/0008094| . 



17 



I. Oda and S. Yahikozawa, "Linking Numbers And Variational Method, 
Phys. Lett. B 238, 272 (1990). 



R. R. Landi m and C. A. S. Almeida, Phys. Lett. B 504, 147 (2001) [hep-| 
|th/0010050| . 

R. R. Landim, Phys. Lett. B 542, 160 (2002) |hep-th/0108242| . 

M. Nakahara, "Geometry, topology and physics," Boca Raton, USA: Taylor 
& Francis (2003) 573 p 

S. Kachru, M. B. Schulz and S. Trivedi, "Moduli stabilization from fluxes in 
a simple IIB orientifold," JHEP 0310, 007 (2003) [arXiv:hep-th/0201028] . 

I. Antoniadis and T. Maillard, "Moduli stabilization from magnetic fluxes in 
type I string theory," Nucl. Phys. B 716, 3 (2005) l [arXiv:hep-th/0412008| . 

V. Balasubramanian, P. Berglund, J. P. Conlon and F. Quevedo, "System- 
atics of Moduli Stabilisation in Calabi-Yau Flux Compactifications," JHEP 
0503, 007 (2005) |arXiv:hep-th/0502058| . 

G. Alencar, R. R. Landim, M. O. Tahim, K. C. Mendes, R. R. Landim, 
M. O. Tahim, R. N. C. F ilho and K. C. Mendes, Europhys. Lett. 93, 10003 
(2011) |arXiv:1009.l"T83| [hep-th]]. 

B. Mukhopadhyaya, S. Sen, S. Sen and S. SenGupta, "Bulk Kalb-Ramond 
field in Randall Sund rum scenario," Phys. Rev. D 70, 066009 (2004) 
|arXiv:hep-th/0403098| 

B. Mukhopadhyaya, S. Sen and S. SenGupta, "Bulk antisymmetric ten- 
sor fields in a Ra ndall-Sundrum model," Phys. Rev. D 76, 121501 (2007) 
|arXiv:0709.3"428| [hep-th]]. 

G. De Risi, "Bouncing cosmology from Kalb-Ramond Braneworld," Phys. 
Rev. D 77, 044030 (2008) [arXiv: 071 1.3781] [hep-th]]. 

B. Mukhopadhyaya, S. Sen and S. SenGupta, "A Randall-Sundrum sce- 
nario with bulk dilaton and torsion," Phys. Rev. D 79. 124029 (2009) 
[arXiv:0903.6722| [hep-th] ] . 

G. Alencar, M. O. Tahim, R. R. Land im, C. R. Muniz and R. N. Costa 
Filho, Phys. Rev. D 82, 104053 (2010) |arXiv:1005.l"691~1 [hep-th]]. 

A. Kehagia s and K. Tamvakis, Phys. Lett. B 504 (2001) 38 [hep-| 
|th/0010112] . 

G. Alencar, R. R. Landim, M. O. T ahim, C. R. Mun iz and R. N. Costa 
Filho, Phys. Lett. B 693, 503 (2010) |arXiv:1008.0678| [hep-th]]. 



R. R. Landim, G. Alencar, M. O. Tahim, M. A. M. Gomes and R. N. Costa 



Filho, Europhys. Lett. 97, 20003 (2012) arXiv:1010.1548 [hep-th]] 



18 



C. -EFu, Y. -X . Liu, K. Yang and S. -W. Wei, JHEP 1210, 060 (2012) 
|arXiv:1207.3152| [hep-th]]. 

N. Arkani-Hamed, S. Dimopoulos, G. R. Dvali and N. Kaloper, JHEP 
0012, 010 (2000) [hep-ph/9911386). 



19 




20 



